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ABSTRACT 

We report on a series of twenty ~10 5 c/s, 0.125 msec time-resolution RXTE obser- 
vations of the Z source and low- mass X-ray binary Scorpius X-l. Twin kilohertz quasi- 
periodic oscillation (QPO) peaks are obvious in nearly all observations. We find that the 
peak separation is not constant, as expected in some beat-frequency models, but instead 
varies from ~310 to ~230 Hz when the centroid frequency of the higher-frequency peak 
varies from ~875 to ~1085Hz. We detect none of the additional QPO peaks at higher 
frequencies predicted in the photon bubble model (PBM), with best-case upper limits 
on the peaks' power ratio of 0.025. We do detect, simultaneously with the kHz QPO, 
additional QPO peaks near 45 and 90 Hz whose frequency increases with mass accretion 
rate. We interpret these as first and second harmonics of the so-called horizontal-branch 
oscillations well known from other Z sources and usually interpreted in terms of the mag- 
netospheric beat-frequency model (BFM). We conclude that the magnetospheric BFM 
and the PBM are now unlikely to explain the kHz QPO in ScoX-1. In order to succeed 
in doing so, any BFM involving the neutron star spin (unseen in ScoX-1) will have to 
postulate at least one additional unseen frequency, beating with the spin to produce one 
of the kHz peaks. 

Subject headings: stars: individual (Sco X-l) — stars: neutron — pulsars: general 
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1. Introduction 

Kilohertz quasi-periodic oscillations (QPO) have 
now been observed in eleven low-mass X-ray binaries, 
the Z sources ScoX-1 (van der Klis et al. 1996a; here- 
after Paper 1), GX 5-1 (van der Klis et al. 1996b) and 
GX17+2 (van der Klis et al. 1997), the atoll sources 
(see Hasinger and van der Klis 1989) 4U 1728-34 
(Strohmayer et al. 1996a), 4U 1608-52 (Berger et al. 
1996), 4U 1636-53 (Zhang et al. 1996), 4U 0614+09 
(Ford et al. 1997), 4U 1735-44 (Wijnands et al. 
1996a) and 4U 1820-30 (Smalc, Zhang and White 
1996), and in KS 1731-260 (Morgan and Smith 1996) 
and a source near the galactic center, perhaps 
MXB 1743-29 (Strohmayer et al. 1996b). In most of 
these sources, the QPO frequency has been observed 
to increase with accretion rate M; frequencies are in 
the range 325-1 193 Hz and relative peak widths vary 
between 0.11% and 10%. Most often double peaks are 
observed, with a separation in the range 250-360 Hz. 
In 4U 1728-34 (Strohmayer et al. 1996a), during 
X-ray bursts, a third peak is seen near a frequency 
of 360 Hz, compatible with the separation of the two 
kHz peaks, which remains constant as the peaks move 
up and down in frequency. Three other cases of 
three similarly commensurate frequencies have been 
reported (4U 0614+09; Ford et al. 1997, 4U 1636-53; 
Zhang et al. 1997, and KS 1731-260; Wijnands and 
van der Klis 1997). This strongly suggests a beat- 
frequency interpretation, with the ~360Hz peak in 
4U 1728—34 at the neutron star spin frequency, the 
higher- frequency (hereafter "upper" ) kHz peak at the 
Kepler frequency corresponding to some preferred or- 
bital radius around the neutron star, and the lower- 
frequency (hereafter "lower") kHz peak at the dif- 
ference frequency between these two. Strohmayer et 
al. (1996a) suggested that this preferred radius is 
the magnetospheric radius. Miller, Lamb and Psaltis 
(1996) proposed it is the sonic radius. 

Although ways out can always be found, this class 
of models naturally predicts the peak separation to 
be constant. In this paper we present data that show 
conclusively that in Sco X-l the peak separation varies 
systematically. A brief announcement of this result 
already appeared in van der Klis et al. (1996c). We 
also present evidence for the presence of horizontal 
branch oscillations (HBO; see van der Klis 1995 for a 
recent summary of Z-source characteristics) in Sco X- 
1 near 45 Hz with a harmonic near 90 Hz. This is the 
first time that HBO have been positively identified 



in Sco X-l. HBO arc usually interpreted in terms of 
the magnetospheric beat frequency model (Alpar and 
Shaham 1985, Lamb et al. 1985), precluding the ap- 
plication of this model to the kHz QPO that occur at 
the same time. 

2. Observations and analysis 

We observed Sco X-l with the RXTE PCA (Bradt, 
Rothschild and Swank 1993) 20 times during May 24- 
28, 1996. Each observation consisted of 2-4 continu- 
ous data intervals of 1-3 ks each. Single- and double- 
event data (Paper 1) were recorded in parallel and 
combined off-line to enhance sensitivity. A time reso- 
lution of 1/8192 s (^0. 125 ms) was used throughout. 

During these observations various offset angles were 
used, and not all 5 detectors were always on. For these 
reasons, the expected Z track in the X-ray color-color 
diagram can not now be recovered; this awaits bet- 
ter understanding of the spectral calibration of the 
PCA at high count rates and off-axis source positions. 
Raw count rates varied between 60 and 130kc/s (2- 
60keV). 
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Fig. 1. — Power spectrum from ~+0 ks of data show- 
ing double kHz QPO peaks, with best fit superim- 
posed. Note the absence of additional peaks. The 
sloping continuum above 1kHz is instrumental (§2). 

We calculated power spectra of all 0.125 ms data 
using 16-s data segments and calculated one average 
spectrum for each continuous data interval. For mea- 
suring the properties of the kHz QPO we fitted the 
256-4096 Hz power spectra (Fig. 1) with a function 
consisting of a constant, two Lorentzian peaks and 
either a broad sine or a broad sinusoid to represent 
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the deadtime-modified Poisson noise, depending on 
the Very Large Event window setting (Zhang et al. 
1995, Zhang 1995). The PCA deadtime process at 
10 5 c/s is not, as yet, sufficiently well understood to 
accurately predict this Poisson component. We can 
not, therefore, report on the properties of any intrin- 
sic broad noise components in the kHz range. 

For measuring the 45 Hz QPO and its harmonic 
we fitted the 16-256 Hz power spectra with a broad 
Lorentzian centered near zero frequency to represent 
the continuum, and one or two Lorentzian peaks to 
model the QPO. The conversion of the power in the 
QPO peaks to fractional rms amplitude depends on 
the derivative of the deadtime transmission function 
with respect to count rate (van der Klis 1989), which 
we do not know. The deadtime is expected to sup- 
press the QPO amplitude more than the total count 
rate. Our reported raw (i.e., uncorrected for dead- 
time) fractional rms amplitudes are therefore lower 
limits to the true values. These could be several times 
as large. 

3. Results 

Kilohcrtz QPO were detected in all observations. 
The peaks (Fig. 1) are very significant, with raw rms 
values of up to 2.5%, and the spectra are well-fit by 
the fit function described in §2. Fig. 2 illustrates the 
changes in power-spectral shape as a function of in- 
ferred M. Notice the increase in frequency and de- 
crease in power of the two kHz QPO peaks, the emer- 
gence of the normal-branch oscillations (NBO) near 
6 Hz apparently from the low- frequency noise (LFN), 
and the complicated variations in strength and shape 
of the 45 and 90 Hz peaks with M (increasing up- 
wards). As in Paper 1, the frequency of the NBO is 
correlated to that of the kHz QPO. 

As we can not estimate M from the X-ray color- 
color diagram, we plot in Fig. 3 the results of our 
fits vs. the centroid frequency v u of the upper peak. 
v u increases monotonically with position along the Z 
track, and, by inference, M (Paper 1). Inspection of 
the variations in the properties of the NBO, whose 
relation to M is known, confirms this. 

Fig. 3 a shows the variation of the kHz peak sep- 
aration. There is a strong decrease with v u : when 
their frequencies increase, the two peaks systemati- 
cally move closer together. The relation of peak sep- 
aration to v u appears to be non-linear. Fig. 36 shows 
that the higher-frequency QPO become more coher- 
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Fig. 2. — Representative ~3ks power spectra sorted 
according to inferred M (increasing upwards) and 
shifted up by factors of 1, 2, 5, 10, 20, 40 and 80, 
respectively, for clarity. The frequency of the upper 
kHz peak increases with M from 872 to 1115 Hz, that 
of the lower one from 565 to 890 Hz. The large width 
of the 10-Hz peak in the top trace is due to peak 
motion. The sloping continua in the kHz range are 
instrumental (§2). 

ent as their frequency increases. The same trend 
is present, less pronounced, in the lower- frequency 
QPO. The ratio of upper to lower kHz-peak width 
(not shown) is roughly constant at 1. Fig. 3 c shows 
how the raw fractional rms amplitude of the upper 
peak falls with M. Again, the lower peak shows a sim- 
ilar trend, but less pronounced. Although the raw rms 
values are deadtime- and therefore count-rate depen- 
dent (§2), the changes in raw count rate are random 
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Fig. 3. — Kilohertz peak separation (a), relative 
width (full width at half maximum divided by cen- 
troid frequency, expressed in percent; b) and raw 
fractional rms amplitude (c) of the highest frequency 
("upper") peak, ratio of the power in the upper to 
that in the lower peak (d), raw count rate (e) and 
horizontal branch oscillation first and second harmon- 
ics' frequency (f) as a function of the frequency of the 
upper kHz peak. One data point at 1075 Hz with a 
value near 3 is off-scale in (d). Eight low signal-to- 
noise data points were not plotted. The gap between 
1000 and 1075 Hz is due to a lack of good data. 

(Fig. 3e), so that the trend in rms must be intrinsic to 
the source. Fig. 3d shows the kHz peaks' power ratio 
(which is essentially free of deadtime effects, as any 
time scale dependence in these effects is expected to 
become appreciable only for time scales near the dead- 
time, — 10/zsec; Zhang 1995). The power ratio varies 
non-monotonically as a function of v u , dropping by 
more than a factor of 2 between 875 and 1000 Hz, 
then increasing again. 

We detect no other kHz QPO peaks. When v u 
<lkHz and the detected peaks are relatively strong, 



any other kHz peaks of similar width are typically less 
than 0.025-0.2 times the power in the detected peaks 
(95% confidence), depending on data selection. For 
^ tt >lkHz these limits worsen as the detected peaks 
get weaker. 

Fig. 3/, finally, shows the frequencies of the ~45 
and ~90Hz peaks. For ^„<960Hz the 45 Hz peak 
has a width of -8 Hz and the 90 Hz peak of 20-40 Hz. 
Both peaks drop rapidly in raw rms amplitude as v u 
increases, from ~1% at j/ u — 850 to below the detection 
limit of -0.4% between 960 and 1000 Hz. There are 
no good data between 1000 and 1075 Hz. At 1075 Hz 
the —45 Hz peak reappears at 1% rms, now much 
broader (20-40 Hz; see also Fig. 2). 

4. Discussion 

Our data place severe constraints on any model 
for kHz QPO so far discussed. We note, that the twin 
kHz QPO in ScoX-1 and those in the atoll sources are 
likely to be the same phenomenon. The frequencies, 
their dependence on M, the coherencies, the peak sep- 
arations and the fact that there are two peaks, one 
of which sometimes becomes undetectable at extreme 
M, are too similar to be attributed to just coinci- 
dence. Even the amplitude of the kHz QPO in Sco X-l 
is in the range of that seen in the atoll sources (al- 
though in some of them much higher values have been 
seen). This implies, then, that the variable peak sep- 
aration we detect in Sco X-l must be explained within 
the same model as the properties of the twin peaks in 
atoll sources. 

The 45 and 90 Hz QPO we reported here are nearly 
certainly the first and second harmonic of the hori- 
zontal branch QPO (HBO), usually seen in Z sources 
in the horizontal and upper normal branches (see van 
dcr Klis 1995). This was already tentatively suggested 
for the broad 45 Hz peak seen when j/ u >1075 Hz in Pa- 
per 1. The similarities with HBO in other Z sources 
include the frequency, its increase (but see Wijnands 
et al. 1996), and the decrease in rms, with M, and 
the presence of a 2nd harmonic. This constitutes the 
first positive identification of HBO in Sco X-l. 

The morphology of the spectra in Fig. 2 seems to 
suggest that the well-known, slightly peaked broad 
noise component below 20 Hz usually called 
low-frequency noise (LFN) "peaks up" into the, also 
well-known, 6 Hz normal branch QPO (NBO) when 
M increases. Appearances may deceive. One way to 
check whether this suggested relation between NBO 
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and LFN is real would be to study the photon en- 
ergy dependence of the amplitude and phase of the 
variability, which in NBO can be quite characteristic 
(Mitsuda and Dotani 1989). There is a small shoulder 
in our spectra below the 45 Hz peak (see, e.g., Fig. 2, 
second spectrum from below) that can perhaps be 
identified as the true signature of the noise compo- 
nent that is expected to accompany HBO (cf. Lamb 
ct al. 1985). 

The photon bubble model (Klein et al. 1996), and 
also some neutron star vibration models, predict sev- 
eral kHz QPO peaks at frequencies above those of the 
two detected ones, and of similar strength as these. 
The fact that we observe just these two kHz peaks, 
with good upper limits on any additional ones, is a 
strong argument against these models. 

If the kHz QPO are due to a millisecond X-ray 
pulsar whose pulsations we see (Doppler-shifted) re- 
flected off inhomogeneities in the Fortner et al. (1989) 
radial flow (Paper 1), then one would expect the QPO 
to become weaker at low M, when the radial flow sub- 
sides; instead we find that the kHz QPO become much 
stronger at low M. In a variant on this model, sup- 
pose that two relativistic jets are emerging in opposite 
directions from near the neutron star, and that we are 
seeing the signal from a central millisecond X-ray pul- 
sar, not directly, but reflected off inhomogeneities in 
the jets. Because there are two jets, there are two 
QPO peaks, at v T = f pu ; se (l — ^/ c )/(l ± ( v / c ) cos9), 
where v is the jets' speed, 9 their angle with the line 
of sight, and v pu i S e the unseen pulse frequency (van 
der Klis 1996). Such a model fits very well to the non- 
linear relation plotted in Fig. 3a, with j/ pu i se =1370Hz 
(which could be twice the spin rate of the neutron 
star) and 0=61°, if v/c decreases from 0.48 to 0.26 
with increasing M. The model predicts that the X- 
ray spectra from the two QPO peaks should show 
similarly different redshifts as the QPO frequencies. 
However, this kinematic model has no way to account 
for the constant peak separations over a large range 
in QPO frequency reported in atoll sources. 

In view of the observations of three commensu- 
rate frequencies in several atoll sources (§1), beat- 
frequency models (BFMs) are the mechanism of choice 
for explaining kHz QPO. The fact that in ScoX-1 
(this paper), GX5— 1 (van der Klis et al. 1996b) and 
GX17+2 (van der Klis ct al. 1997) we observe HBO 
and kHz QPO simultaneously shows conclusively that 
these can not both be explained by the magnetospheric 
BFM. If this mechanism produces the HBO (§1), then 



the kHz QPO need another model. 

Beat-frequency models with the neutron star spin 
as one of the participating frequencies predict a con- 
stant kHz-peak separation, which so far was consis- 
tent with observations. However, this prediction is 
clearly contradicted by our ScoX-1 data. In order to 
explain the data, such BFMs would have to be modi- 
fied such that, in addition to the unseen (in ScoX-1) 
spin frequency v s there is another unseen frequency 
beating with v s to produce one of the two kHz peaks. 
The sonic-point beat-frequency model (Miller, Lamb 
and Psaltis 1996) may allow such modification (Lamb 
1996, priv. comm.). If this model explains the kHz 
QPO and the magnetospheric beat-frequency model 
the HBO, then in ScoX-1 the sonic radius is approx- 
imately [v u /{vhbo + i / s )] _2 ^ 3 ^ j 0.5 times the magne- 
tospheric radius, implying the presence of a consid- 
erable near-Keplerian flow, where clumps remain in 
stable orbit for up to 10 2 cycles, well within the mag- 
netosphere. 
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